The race for supersymmetry: using m T2 for discovery 
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We describe how one may employ a very simple event selection, using only the kinematic variable 
77iT2, to search for new particles at the LHC. The method is useful when searching for evidence 
of models (such as i?-parity conserving supersymmetry) which have a Zi parity and a weakly- 
interacting lightest parity-odd particle. We discuss the kinematic properties which make this variable 
an excellent discriminant against the great majority of Standard Model backgrounds. Monte Carlo 
simulations suggest that this approach could be used to discover supersymmetry with somewhat 
smaller integrated luminosities (or perhaps lower center-of-mass energies) than would be required 
for other comparable analyses. 



I. INTRODUCTION 

The search for new particles cannot be divorced from 
measurement of their masses. Indeed new particles have 
usually been discovered by measuring differential distri- 
butions in variables sensitive to their mass [H, ■ A no- 
table example was the successful W boson search at the 
CERNSppS collider using the transverse mass, tot- 
Using the right variables can improve the discrimination 
between signal and background, increase the statistical 
power and thus expedite discovery. 

One of the primary goals of the Large Hadron Col- 
lider (LHC) is to search for physics beyond the Stan- 
dard Model. Many well-motivated models of new physics 
0| have a Z2 parity (for example i?-parity in supersym- 
metry) which predicts the presence of stable, presumed 
weakly-interacting, particles which would contribute to 
the astronomical dark matter. Such particles must be 
produced in pairs and would be invisible to the LHC de- 
tectors. 

In this paper we explain why the 777,^2 variable, which 
was originally proposed [5( for supersymmetric particle 
mass determination in hadron colliders, is also useful in 
searches for such models. 1 

We start by recalling the definition of the transverse 
mass and its generalisation tot2 j and outline their defin- 
ing properties (Section [H]) . In Section UTT1 we show that 
experimentally-interesting kinematic configurations of- 
ten lead to near-minimal values of tot2- In Section ITVl we 
explain why these properties should cleanly separate new 
physics signals from the great majority of the Standard 
Model background events. A Monte Carlo simulation of 
a typical di-jet channel can be found in Section [V] We 
discuss the results and conclude in Section fVTl 
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1 Some of these observations were originally made in a talk at 
CERN by the authors in May 2007. Since then, the method has 
been employed, for example in [f|, but the underlying proper- 
ties which represented the motivation for the approach have not 
previously been published. 



II. DEFINITIONS 

The familiar transverse mass, which was used for the 
W boson discovery, is defined by 



to t 
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where m v and rzij are the masses of the visible and invis- 
ible systems respectively (e.g. the electron and neutrino 
for the W case), vj- and qr are their momenta in the 
transverse plane, and the transverse energies are defined 
by 



and 



The transverse momentum of the unobserved neutrino 
can be inferred from momentum conservation in the di- 
rections perpendicular to the beam: it is assumed to be 
equal to the missing transverse momentum, p which is 
defined to be the the negative sum of the transverse mo- 
menta of all the particles observed in the detector (i.e. 
including hadrons in the W case). 

Why was mj- useful for the W search? The defining 
property of Tot is as follows: for events in which a visible 
and an invisible system originate from the decay of a 
single parent particle of mass mo, and when the correct 
daughter masses are used, mx satisfies 



rriT < Too 



(2) 



by construction 2 , with equality when the rapidity of all 
invisible daughters is equal to that of the visible system. 
The transverse mass gives the lowest event- by- event up- 
per bound on tuq. The maximum value of Tot over all 
events therefore delineates the boundary between allowed 
and forbidden values of the parent mass, mo- 

For the W boson search the allowed values of tot could 
therefore span a range m< < tot < mw for signal 
events. The lower bound is given by m< = m, + to,-, 
which is zero if one assumes (as was done implicitly in 



2 If the width of the particle and the experimental resolutions are 
assumed to be small. 
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the original W search papers) that both m v and rrii are 
zero. Other processes which could have led to lepton + 
p T final states, for example leptonic decays of r leptons 
or B hadrons were also constrained by similar inequal- 
ities, but these backgrounds have smaller values of mo- 
There was therefore a region bounded approximately by 
m B < mT < "iff in which there were very few back- 
ground events and a preponderance of the signal. 3 

A corresponding construction can be made for the 
types of LHC search described in the introduction. We 
are now interested in the case in which two new particles 
are produced, each of which decays to a set of (one or 
more) visible and (one or more) invisible daughters. We 
label the two branches of the decay with the superscripts 
W and ( 2 ) to distinguish them, and for the purposes of 
this paper assume that the visible particles can be unam- 
biguously assigned to one or other parent, though this 
need not be done in general Q. According to our as- 
sumptions, each invisible system should include one dark 
matter candidate. 

As was the case for single-particle decays, there will be 
some values of parent and (visible and invisible) daughter 
masses which are consistent with the observed momenta, 
and others which are not. 

The p T constraint needs some modification compared 
to the single decay case. For a pair of invisible sys- 
tems the transverse momentum of each the two invisible- 
particle systems is not individually known, but the sum 

is constrained by + q^ = p T . One can construct 
the boundary between the allowed and disallowed regions 
using the variable tot2 [1] 

/ „, (1) (2)s f ( (1) (2) N 

tot2 (^l i V2 , V> T , to^ ,m\ ) = mm < max I m T , m T 

(3) 

For decays of two parents, each of mass too, to^2 gives 
the lowest upper bound on too consistent with the ob- 
served momenta and the hypothesised daughter masses. 
The maximum of m^i over events therefore delineates 
the boundary of the domain (in the space of parent and 
daughter masses) which is consistent with the observed 
momenta [8]. It is this bounding property which makes 
it ideally suited for measuring sup ersymmetric particle 
masses at the LHC @, H, S El EH ■ It has also recently 
been used to measure the top quark mass at the Tevatron 

a 

We note that the Lorentz 2+1 dimensional vectors 
V = (e^vr) and i = (e^qr) can each represent sin- 
gle particles or multi-body systems. For example v could 
represent an individual lepton, or a jet or a di-jet pair. 
Similarly each of the invisible-particle systems may rep- 
resent a single massive invisible particle (such as ax?) or 
a single particle of negligible mass (neutrino) or a set of 



3 In Q the lower end of this distribution was removed because of 
threshold cuts on the electron px and p. 



invisible particles (e.g. \i + v )- F° r composite systems, 
the composite 2+1 vectors can be formed from the vector 
sum of the 2+1 vectors of their individual constituents. 4 

III. RELEVANT PROPERTIES 

Later in Section IIVI and [Vj we shall find that many 
background processes predict values of tot2 near the 
global minimum, m<. The value of this minimum is 
therefore of some interest. It follows from ([3]) that 

to< = max(m^',m^) , (4) 

where m< = ml + m\ is the global minimum of each 

of the two individual . 

In order to construct tot2 we require hypotheses for the 
invisible particle(s) masses m^ n \ In a search we usually 
don't know a priori whether any invisible particles are 
massive (or if so, what those masses are) so we choose 
to set the invisible-particle mass hypotheses to.; to zero. 
This is the only value guaranteed to preserve © if the 
invisible particles' true masses are unknown. 5 It is also 
the appropriate choice for background processes where 
the invisible particles are (~ massless) neutrinos. 

In the rest of this section we show that the observ- 
able tot2(^i, «2, 0, 0) must adopt small values for a 
variety of kinematic configurations. In Section IIVI we 
use these results to show that such configurations are 
satisfied (at least approximately) by very many of the 
Standard Model processes which represent backgrounds 
to example searches. 

Lemma 1 When a pair of particles are produced, both 
with mass too, and each parent decays to a visible system 
v and an invisible system i, then tot2(^i, «2, ]^ T , 0, 0) < 
too ■ 

Proof Consider the chained pair of inequalities 
m T2 (vi,V2,p T ,0, 0) < m T 2(vi,V2,f T ,mf\mf ) ) < m . 
The first inequality follows since each mi^ is a mono- 
tonic function of the (non-negative) parameter m\ . 
The second inequality is satisfied by tot2 by construc- 
tion: one of the partitions of the missing momentum is 
the correct one, and for that partition each tot < too by 

©• 

Lemma 2 When two particles are produced with dif- 
ferent masses m\ and mi and each parent decays to 
a visible system v and an invisible system i then 
tot2(wi,«2,^ t ,0, 0) < max(mi,TO 2 ). 



4 Alternative combinatorial procedures are available for the visible 
system [T3| . 

5 If and when signals suggesting new invisible particles are ob- 
served, their masses can later be determined, for example by 
constructing variables sensitive to the kink in the maximum of 

is as a function of the input value of rrii lldll . 
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Proof As for Lemma[T]except that at the correct p T par- 
tition we can only be sure that mlp < mi and rrSp < mi 

and so max(m^ , ) < max(mi, 7712) for that parti- 
tion. 

Lemma 3 When ' = and rriv — m\ — then 
iriT2 — m<; for n E {1, 2}. 

Proof Without loss in generality let n = 2 (and so n = 

(2) (2) (2) 

1). Now m T = m< V . There exists a partition 
of Tp T with c£p = for which vrSp — m< . The result 
follows from © and ©. 

An important example of Lemma [3J is the decay of 
a single parent of mass mo. The second (non-existent) 
visible system is represented by V2 = (0, 0). With m\ — 

0, Lemma [3J shows that m T2 = m< = max(mf \ mi 1} ) 
where the second equality is a result of (|4|). 

(1 2) 

Lemma 4 When ip T = and mj = i/ien mT2 = 
m< . 

Proof For p T = there exists a trivial partition of the 
missing momentum with q^! ' — qS' = 0. For that par- 
tition, rrSp = m<^ and = m^; the result follows 
from (J3J) and ©. 

Lemma 5 When m^' 2 ^ = 0, mi™' = and || 
then m,T2 = m< for n £ {1,2}. 

Proof Without loss of generality let n = 1. There exists 
a partition of ip T with q^' = Tp T and q^ = 0. Each 
tut is equal to its global minimum by |T]) , and the result 
follows from 

Lemma 6 When mj 1 ' 2 '' = m« = 0, and |A can be 

expressed as the sum p T = x^qj^ for some real non- 
negative Xk then vriT2 = rn < = 0. 

Proof For the partition of jA given by these Xk, q^ || 
Vy simultaneously for both n 6 {1,2}. For that parti- 
tion each = m<^ = 0. The result follows from Q 
and (gj). 

Lemma 7 J/ either or both visible systems are composite, 
then each of the results in Lemmas [7] - also hold when 
either or both the composite Lorentz 2+1 vectors v are 
replaced by any subset of their ( respective ) constituents. 

Proof Lemmas [3J - [5] follow exactly as before. For Lem- 
mas [1] and [5] it is sufficient to show that tot2 cannot be 
increased when a constituent is removed from a compos- 
ite (visible) system. 

Now 77i|, is the inner product of a sum of Lorentz 2+1 
vectors 

m| = (v + if = c fe) ' 



where the sum runs over all the constitutent Ck , whether 
visible or invisible. Separating one of the visible particles, 
Cj and letting a = J2k^j c fc> 

m\ — a 2 + Vj + 2(T • Vj . 

Since each of the constituent vectors is time-like (or null) 
each of the terms is non-negative and a 2 < m\. But 
Va^ is precisely the transverse mass one obtains if Vj is 
omitted, so the transverse mass cannot be increased by 
omitting a particle. This result also holds at the parti- 
tion chosen by the minimisation of (J3J so m,T2 is never 
increased by omitting one (or by induction more than 
one) of the visible particles. 

Physically relevant configurations will rarely (if ever) 
conform to the precise configurations described in Lem- 
mas [U - [71 parent particles will be off-mass-shell, Stan- 
dard Model particles will have small but non-zero masses, 
p will never be exactly zero nor vectors exactly parallel. 
However, since is a continuous function of its inputs, 
kinematic configurations close to those described above 
will have upper bounds close to these idealised cases. 

IV. EXAMPLE 

The properties described in Section |TTT] turn out to be 
particularly useful for various LHC searches. We demon- 
strate the reasons why by examining the concrete ex- 
ample of pair-production of heavy, strongly-interacting 
particles of mass mo, each of which decays to a light- 
quark jet and an invisible particle. The characteristic 
final state is thus two (usually high-p^) jets with signifi- 
cant p . Examples of models which could lead to such a 
final state can be found in lTablcTl 

Supersymmetry : q q — > q\i qXi 

UED : qi qi — > 071 571 

Leptoquarks : LQ LQ — > qvqv. 

TABLE I: Examples of models of new physics producing di- 
jets in association with p . 

When considering how to search in such channels, a 
variety of levels of sophistication can be envisioned. 

In a typical cut-based search one would select events 
with large p T and two high-p<r jets; require sufficiently 
large azimuthal angles A(f)j between each jet and the p 
vector to reduce backgrounds from "fake" (from the mis- 
measurement of one of the jet energies) or "real" (from 
neutrinos in jets) p from the jet; and apply further cuts to 
reduce background processes (such as leptonic ti decays) 
which produce neutrinos and high-p-p jets with larger 
Atpj Q . The remaining events should be signal-enriched, 
and so, provided that the residual background contribu- 
tions can be understood, one can check for an excess 
of events not explained by the Standard Model. The 
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Process 


m T 2(«i,U2, ^ T ,0, 0) 


Comments 


QCD di-jct — > hadrons 


= max mj by Lemmas 1 1 141 




QCD multi jets — ► hadrons 


= max rrij by Lemma [4] 




ti production 


= max m.j by Lemma [4] 


fully hadronic decays 




< mj by Lemmas 1 1171 


any leptonic decays 


Single top / tW 


= max m j by Lemma [4] 


fully hadronic decays 




< mt by Lemmas 12171 


any leptonic decays 


Multi jets: "fake" f 


= max rrij by Lemma [5] 


single mismeasured jet" 




= max rrij by Lemma [6] 


two mismeasured jets" 


Multi jets: "real" f 


= max mj by Lemma [5] 


single jet with leptonic b decay" 




= max mj by Lemma [6] 


two jets with leptonic b decays" 


Z — > vu 


= by Lemma [3] 




Zj -> uuj 


= mj by Lemma [3] 


one ISR jet" 


W tv b 


= by Lemma [3] 




W j-> £uj b 


< mw by Lemma [2] 


one ISR jet" 


WW -> iviv b 


< mw by Lemma [1] 




ZZ — ■+ vvvv 


= by Lemma El 


also = rrij for one ISR jet" 


LQ LQ — * qvqv 


< rriLQ 




qq^ qxi<ix° 


< rriq 


^" i.e. can take large values 


qi,qi -> 97i>97i 


< m qi 





° Assuming that the relevant jet(s) are identified with one of the two visible particle systems v^ n \ 
b Even if the lepton is mistaken for a jet. 



TABLE II: Examples of Standard Model backgrounds, various signals, and associated values of 
m,T2- In the cases marked 6 , there are more restrictive upper bounds if the lepton is assumed 
not to be misidentified as a jet. 



most difficult task is to understand the contribution of 
the background to the selected events - for which a var- 
ious techniques have been proposed [f| [lij]- Applying 
such a set of many single- variable cuts is fairly straight- 
forward but is somewhat wasteful of signal events (which 
one can assume will be in scarce supply at the time of 
any interesting discovery). 

At the other extreme, the formally optimal search 
method (needing the fewest events) would require cal- 
culation of the likelihood for all events for every signal 
and background hypothesis. This is often prohibitively 
difficult in practice even if all such hypotheses can be 
identified and modelled. 

A pragmatic approach, intermediate in complexity, is 
to find a single, easily-calculable observable which, while 
not optimal in the formal statistical sense, still gives very 
good discrimination between the majority of signal-like 
and background-like events, based on some rather general 
principles such as relativistic kinematics. In what follows 
we shall find that tot2 is just such an observable. 

For the example search channels, we identify the two 
visible systems v^ 1 ' with the two highest-px jets. Our 
jets, though massive, should have a sufficiently small 
mass that rrij ~ is a good approximation. For the 

reasons discussed in Section Hill we set both nS™' to zero. 
For these choices, m< is small (equal to the larger of the 
) and the mass conditions required for Lemmas [3] - 
[7] are satisfied in the rrij — > limit. 



ITablc III lists a variety of Standard Model processes 
which form backgrounds to the search channel. From 
Lemmas [T]- [7] we expect to find restrictive upper bounds 
on rriTi for many of these. Processes with small p T ; 
with small jet pr; with small A0j; or from production 
of one or two Standard Model particles; all bound rriTi 
from above. Such processes constitute the majority of 
both the "physics" and "detector" backgrounds to this 
channel so one can remove the vast majority of the back- 
ground simply by requiring large rriT2- No additional 
cuts are required - other than perhaps modest trigger 
requirements on jet pt or on p . 

An example "physics" background is ti production. 
This is pair-production of equal-mass particles, so rriT2 
is bounded above by the mass of the top quark to* by 
Lemma [TJ One could argue that this is not a very strict 
bound, but it is much smaller than masses of the par- 
ticles typically being searched for, and the top is the 
heaviest known Standard Model particle, so other sim- 
ilar background processes, e.g. hadronic t + t~ decays 
should adopt smaller values of ttiti again. 

An example of a detector-induced background is "fake" 
p T - events with no invisible particles, but where a sub- 
stantial fraction of the energy of one of the leading jets is 
lost. This can happen when energy is deposited in inac- 
tive material the detector (e.g. cables, services or support 
structures). Such pathological events cannot necessarily 
be cleanly identified. In these cases the detector usu- 
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ally measures a jet with Lorentz 2+1 vector j' = aj 
(0 < a < 1), and one gains a contribution to ]f> T of 
(a — 1)jt- In the absence of any other source of missing 
momentum, || jr, so mT2 — * to< by Lemma [5] or |6] 
Similar arguments apply to heavy-quark jets where lep- 
tonic decays lead to production of neutrinos close to the 
jet axis. 

The other backgrounds in I Tabic III are also forced to 
small values of to*t 2 - The least restrictive is < m t since 
the top is (assumed here to be) the heaviest Standard 
Model particle. 

What values of tot2 are expected for any new parti- 
cles? At the upper end it is clear from Lemma [TJ that 
tot2 _ m-0 for the processes in lTablc HI One needs a sig- 
nificant number of events with rriT2 > m t to have a signal 
region which is relatively free of background. Now if the 
correct value of rrii were to be used then the upper bound 
(mo) would be saturated since there is a significant den- 
sity of states with rriT2(vi, i>2, p T , mi, mi) close to mo. 
We have chosen to input the lowest conceivable value 
rrii — > rather than the true invisible-particle mass, so 
the argument does not prove that the bound is saturated. 
However one can see from ([1} that, provided rrii <C \ckr\, 
then events which are close to maximal when the true 
rrii is used will also remain close when we replace this by 
rrii = 0. 

We therefore expect to find a large number of signal 
events, and very little background, in the region approx- 
imately bounded by m t < 771x2 < mo, where mo is the 
mass of the new particle. 



V. SIMULATION 

To illustrate the results of Section IIIII and Section IIVI 
we generate Monte Carlo signal and background samples 
with Herwig++ 2.3.2 (la |. The background processes 
simulated are QCD, it, W — ► £^+jets, Z — > £ + £~+]ets, 
and Z — ► vD+iets. The contribution from diboson+jets 
is expected to be very small @ so is not considered 
here. In order to generate sufficient QCD events in the 
high-pT region, eight samples were generated in slices 
of the px of the hard scatter. For the SUSY signal, 
the SPSla point [l6[ is used (m = 100 GeV, m l / 2 — 
250 GeV, A = - 100 GeV tan /3 =10, n > 0), as cal- 
culated by SPheno 2.2.3 [oj. Table fTTTI lists the lead- 
ing order cross sections calculated by Herwig++, and 
the number of events generated for each of the processes 
considered. 

We cluster hadrons (and 7r°s) with fiducial pseudora- 
pidity (\rj\ < 5) and momentum (j>t > 0.5 GeV) into jets 
using the fastjet [l8| implementation of the anti-^T al- 
gorithm [l|l] . We use the E combination scheme and set 
R = 0.4 and pf in = 10 GeV. To simulate the detector 
effects we smear the majority (1 — e) of the jet energies 




100 200 300 400 500 600 

m T2 [GeV] 



FIG. 1: Distribution of ot,t2 f° r events with two or more jets 
with pr > 50 GeV (and no other cuts). For the signal point 
the squark masses are in the range 500 < rriq < 600 GeV and 
the gluino mass is close to 600 GeV. 



by a Gaussian probability density function of width 

a{E)/Ej = (o.S/^Ej [ GeV] J ® 0.03 

where Ej is the unsmeared jet energy. This resolution 
is typical of one of the general-purpose LHC detectors 
0, [lj|- Since the tails of the distribution are of- 
ten dominated by badly mismeasured jets, we simulate 
pathological energy-loss by applying a different smearing 
function to the remaining fraction (e = 0.1%) of the jets 6 
with probability density: 

p{e) = \ 2E/e2 i {or (°< E < E j) 

[ elsewhere 

The missing transverse momentum is calculated from the 
negative vector sum of the visible fiducial hadrons (in- 
cluding 7T°) and is corrected for the jet smearing. We 
impose the simple requirement that each event contains 
at least two jets with pr > 50 GeV. We then take the two 
highest pt jets as ji t 2 and calculate TOT2O1, ji, ]/> T , 0, 0), 
for all events. We normalise to 100 pb _1 (using the lead- 
ing order cross sections for both signal and background) . 
The resulting distribution can be seen in |Figure l] 



6 @1 suggests a larger value of e ~ 1%. We find a smaller value 
better matches the tails of the p T and mj-2 distributions found 
in full simulation. The detailed form of the transfer function 
clearly needs to be determined from the collision data, but our 
findings are not materially altered by changing epsilon from 0.1% 
to 1%. 
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sample 


LO cross section (pb) 


number of events 


QCD (17 <p T < 35GeV) 


9.44 x 10 s 


6 x 10 7 


QCD (35 < p T < 70 GeV) 


5.99 x 10 7 


5 x 10 7 


QCD (70 <p T < 140 GeV) 


3.45 x 10 6 


5 x 10 7 


QCD (140 <p T < 280 GeV) 


1.57 x 10 5 


5 x 10 7 


QCD (280 < PT < 560 GeV) 


5280 


4 x 10 7 


QCD (560 <p T < 1120 GeV) 


116 


3 x 10 7 


QCD (1120 < p T < 2240 GeV) 


1.11 


2 x 10 7 


QCD (p T > 2240 GeV) 


1.15 x 10" 3 


5 x 10 6 


tt 


231 


2 x 10 7 


W -> ^+jets 


17,100 


1.5 x 10 s 


Z -f e+i^+jets 


1780 


9 x 10 7 


Z — » vv+jets 


3440 


5 x 10 7 


SUSY signal (SPSla) 


13.5 


2 x 10 6 



TABLE III: Signal and background processes together with their leading order cross-sections, 
and number of events generated for each. 



As predicted, the region m t < < m q is domi- 

nated by signal events. The great majority of the Stan- 
dard Model background events are found at small m^, 
as required by the arguments of Section ITVl The remain- 
ing backgrounds for which tuti is not bounded above 
by the arguments in Section IIVI are: Z — > vv in asso- 
ciation with two hard jets from initial state radiation 
(ISR); W — > leptons plus two hard ISR jets; and three- 
(or more-) jet production where one of the jets looses a 
very large amount of energy, and that jet is not one of the 
two highest px jets input to m^. For the signal point 
examined, and with the level of sophistication used for 
our simulations, we find that the residual backgrounds 
for m,T2 > Trit are well below the supersymmetric signal 
predictions. 



VI. DISCUSSION 

The numerical simulations of Section [V] confirm the 
analytic results of Section HTT1 m^2 adopts small values - 
equal to or less than the mass of Standard Model parti- 
cles - for the vast majority of events from the background 
processes. The region tot2 > mt is signal-dominated; for 
SPSla a single cut requiring tot2 > 230 GeV results in a 
statistical significance S/^/S + B of 15, for 100 pb -1 of 
integrated luminosity at ^/S = 10 TeV. This is a higher 
significance than is found when applying to the same 
simulated events the selections of comparable multi-cut- 
based analyses e.g. d, I20I. [2lT|. 

Applying a single-variable kinematic selection based 
solely on ttiti would be sufficient to separate the signal 
from the background. In practice tot2 is unlikely to be 
calculated by the trigger algorithms (certainly not at the 
first level) so some thresholds for the jet ptS and the f 
will be required in order to keep trigger rates within the 
available bandwidth, particularly at higher instantaneous 
luminosities. Our simulations show that adding typical 



trigger requirements has the effect of removing events 
at small tot2 (as would be expected from Lemmas 13141) . 
Reasonable thresholds leave the large witi region - and 
hence the statistical significance - practically unchanged. 
The lower vuti end of the distribution can be recovered 
using lower-threshold prescaled triggers, so it can be used 
as a control region. 

The remaining backgrounds at large witi are domi- 
nated by states originating from more than two (massive 
or high-py) Standard Model particles (e.g. Zjj, Wjj, 
ttj, jjj). The bounds of Section ILTTl do not apply to these 
since the variable niT2 was developed for the explicit two- 
parent case. To reduce the residual contribution from 
these n > 2-particle topologies one might consider us- 
ing, rather than ttiti-, an n-parent generalisation of the 
transverse mass: 

rriTn = nrin {maxm^ } with i G {1, . . . , n}. (5) 

2~2q T =$ T 

Using §5§ n-body background topologies (and also m < 
n-body topologies by a generalisation of Lemma[2|) would 
be bounded from above by the mass of the heaviest par- 
ent. The problem would be that one would no longer 
expect pair-produced signal topologies (such as those in 
ITablc ip to obtain rriTn values close to the mass of the 
new heavy parent. Indeed a generalisation of Lemma [3] 
shows that for the simplest Z2 signal case, di-jet + |* T , 
rriT3 is forced towards the small value (maxnij) so one 
loses the discrimination power of m,T2- Such n > 2 gener- 
alisations are therefore only likely to be appropriate when 
n heavy signal particles are expected to be produced. 

Our introduction focused on the simplest decay topolo- 
gies of pair-produced heavy particles such as -> 
QXi QXi but our simulations also show excellent discrim- 
ination in more complicated cases using the same tot2 
variable. Decay sequences with many steps, or indeed in- 
dividual multi-daughter decays, such as gg — > qqXiI'lXi 
(via three-body decays), could have been input to tot2 
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in a number of ways. There is no unique choice for 
constructing the two 'visible particle systems' from the 
(k > 2) decay products. We could have chosen to form 
two composite systems (two di-jet systems for the gluino 
example), and use these as the visible inputs to m,T2 
0, 0|- This construction would have provided a large 
number of signal events close to the kinematic bound- 
ary and so would be appropriate for mass determination. 
However forming di-jet composite objects produces visi- 
ble systems which no longer have m, ~ 0, so one would 
lose the desirable properties of Lemmas [3j[5] -0 For our 
search, even though cascade or multi-body decays are ex- 
pected, we have still choosen to form rriT2(ji , 32, p T , 0, 0) 
using only the two highest pt jets — precisely as for the 
simpler di-jet topology. This way we retain the desirable 
properties (Lemmas |31 [5] - [7]) for the backgrounds. The 
m,T2 distribution for the signal can still extend up to large 
values (close to rrig for the three-body decay) , albeit with 
fewer near-maximal events than would be found by using 
composite two-jet visible systems. By inputting only the 
highest two pt jets to tot2, not only do we take advan- 
tage of the background rejection properties, but we also 
combine many signal channels together, forming a larger 
sample of signal events (which therefore enjoys a larger 
statistical significance) . 

Of course even having found a good discriminating 
variable one is not absolved from the need to understand 
the residual background contributions. Previous studies 
0, HH show that the rates of a wide variety of Standard 
Model processes can be measured in control regions using 

7 It has recently been suggested that ISR jets could themselves 
provide a good indicator of mass scale [22ll , 



the LHC collision data. These are then extrapolated to 
the signal region using techniques which themselves are 
validated by Monte Carlo simulation. Our 777^2 signal re- 
gion has a non-negligible contribution from vector boson 
production in association with two (or more) jets from 
initial state radiation. The dependence of such ISR jets 
on parameters such as y/s and 77 is therefore of consider- 
able interest; jets at large |ry| could, for example, provide 
an ISR-dominated control region. 7 The fact that there 
is a large number of background events in the low 7717^ 
region (< m t ) might actually turn out to be rather help- 
ful - one can measure the total background contributions 
from this control region with high precision. 

Complementary measurements will clearly be needed 
to disentangle the various background processes, so many 
observables will be used by the LHC experiments. Our 
analysis and simulations suggest that 777^2, the natural 
kinematic observable for pair- produced particles, ought 
to be at the front of the queue. 
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